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A METHOD FOR CALCULATING HEAT TRANSFER IN THE LAMINAR FLOW REGION OF 

BODIES 

By H. Julian Allen and Bonne C. Look 



SUMMARY 

This report has been prepared to provide a practical method 
far determining the chordwise distribution of the rate of heat 
transfer from the surface of a wing or body of revolution to air. 
The method is limited in use to the determination of heat transfer 
from the forward section of such bodies when the flow is laminar. 
A comparison of the calculated average heat-transfer coefficient 
for the nose section of the wing of a Lockheed 1S-A airplane 
with that experimentally determined shows a satisfactory agree- 
ment. A sample calculation is appended. 

INTRODUCTION 

With the advent of thermal ice-prevention equipment for 
aircraft, the problem of determining the heat which may be 
transferred from the surface of aerodynamic bodies has 
become of considerable interest. Experimental investiga- 
tions (reference 1) have indicated that adequate heating of 
the forward 10 to 20 percent of an airfoil will prevent ice 
formation on the entire airfoil. Moreover, experience has 
shown that in the cruising flight condition in the absence of 
icing conditions if sufficient heat can be supplied to this 
section of the airfoil to raise its surface temperature from 
70° to 100° F. above the temperature of the ambient air, 
ice will not collect on the airfoil. 

To determine whether this surface temperature rise can 
be obtained in the design of a particular installation, it is 
necessary to dete rmin e the rate of heat transfer from the 
airfoil surface to the atmosphere. In the usual cruising- 
flight condition the nose section of an airfoil experiences 
laminar flow. . In the present report a method for calculating 
the rate of heat transfer from an airfoil section subject to 
laminar flow is developed. The method is made general to 
include laminar flows occurring over surfaces subjected to 
positive or negative pressure gradients and may be applied 
to two-dimensional bodies as well as bodies of revolution. 

THEORY 

The following development of an expression for the heat 
transfer between solids and fluids has been based upon 
Reynolds analogy. This analogy was presented by Reyn- 
olds in an early paper (reference 2) in which he suggested 
that, in a fluid, momentum and heat are transferred in the 
same way, and concluded that in geometrically similar 
systems a simple proportionality relation exists between 
heat transfer and fluid friction. 



The analogy may be applied in the case of the heat transfer 
between solids and fluids when the fluid flow is laminar, 
provided the Prandtl number for the fluid is unity. The 
Prandtl number is the dimensionless parameter 

Pr-!f 

where 

c P specific heat of the fluid at constant pressure 

ix absolute viscosity of the fluid 

k thermal conductivity of the fluid 

For air the value of the Prandtl number is 0.73 rather 
than unity but, as is discussed later in this report, experi- 
mental investigations have shown that heat transfer is only 
slightly affected by variation of the Prandtl number so that 
for most practical purposes the Reynolds analogy may be 
applied. The analogy leads to the equation 



h=-^CtP<hVi 



(1) 



where 
h 



heat transfer coefficient defined as the heat trans- 
ferred per unit time from a unit surface area for 
1° difference of temperature between the surface 
and the fluid "outside" the boundary layer 
p fluid density 

Vi stream velocity just outside the boundary layer 
c f surface friction coefficient defined as the frictional 
drag per unit surface area in terms of the local 
fluid dynamic pressure, <ft 
" The surface friction coefficient may be expressed as 



(2a) 



t 2t 

where r is the frictional force per unit area. 

For laminar flow the frictional force r is the product of the 
absolute viscosity and the velocity gradient in the fluid 
boundary layer at the solid surface. For any given bound- 
ary-layer velocity profile, this velocity gradient is directly 
proportional to the velocity Vi and inversely proportional 
to the boundary-layer thickness. It follows that 



2X/i 



(2b) 



where 5 is the boundary-layer thickness and X is a constant 
dependent on the shape of the boundary-layer velocity 
profile and on the definition of 5. It has become customary 
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(reference 3) to define 5 for laminar boundary layers as the 
distance from the solid surface to a point in the boundary 
layer where the dynamic pressure is one-half its value out- 
side the boundary layer. 

Inserting the value of C/ from equation (2b) into equation 
(1) the heat-transfer coefficient becomes 



*-<¥) 



Rearranging and dividing both sides by the thermal con- 
ductivity of the fluid k this becomes 



(3) 



since the Prandtl number must be assumed to be unity by 
the analogy. The parameter (M)/k is nondhnensional and 
might be properly termed "the boundary-layer Nusselt 
number." 

The appearance of Pr in equation (3) might suggest that 
its value for air (0.73) should be substituted to obtain a 
better approximation of the value of the boundary-layer 
Nusselt number. However, the meager experimental evi- 
dence available (reference 4, p. 249) indicates that for 
laminar flows, the Nusselt number is proportional to (Pr)i 
so that using unity for the value of Pr (equivalent to (Pr)°) 
would appear to be preferable. 

In order to determine h from equation (3) it is necessary to 
know the value of X which is a function only of the shape of 
the boundary-layer velocity profile. Numerous experiments 
with airfoils have shown that the velocity profile of the 
laminar boundary layer in the presence of favorable pressure 
gradient (i. e., where the surface pressure gradient is negative 
proceeding in the downstream direction) is closely approxi- 
mated by the Blasius distribution for the flat plate. For the 
Blasius type velocity profile the value of X is 0.765. 

At points on the surface of a body downstream of the 
minimum-pressure point a laminar boundary layer exhibits 
a tendency to separate. Since the velocity gradient at the 
surface decreases as separation develops the value of X must 
diminish until at the separation point its value is zero. It is 
considered that reducing X linearly from the minimum pres- 
sure point to the separation point will satisfactorily approxi- 
mate the actual case, particularly since this region is of little 
importance from the viewpoint of heat transfer. The location 
of the laminar separation point may be calculated by the 
method of reference 5. 

"With the known variation of X along the surface, the heat- 
transfer coefficient may be calculated when the boundary- 
layer thickness has been determined. The value of 5 may be 
found for the region from the stagnation point to the mini- 
mum pressure point by the method of reference 3. Re- 
arranged in a form more convenient for calculation, for air- 
foils the equation is 



a 5.3c 2 / gi/c \ 



.^TWJ 



(4) 



where 



airfoil chord length 

distance along the surface from the stagnation point 



Si distance to the point for which the boundary- 
layer thickness is being computed 
V 0 free-stream velocity 
V velocity outside the boundary layer at s/c 
Vi velocity outside the boundary layer at sjc 

cVa 

R„ body Reynolds number =— 
For bodies of revolution the corresponding equation is 



5.3XY 8i/L \ 



(5) 



where, in addition 
L body length 
r radius of revolution at s 
ri radius of revolution at «j 



LV 0 

body Reynolds number = 



The value of 5 is indeterminate at the stagnation point. 
It is not necessarily zero as may be shown by the following 
approximate analysis. The initial flow near the stagnation 
point on an airfoil will be approximately that over a circular 
cylinder having a radius equal to the radius of curvature 
at stagnation f. On the surface of a circular cylinder the 
velocity at 6, the angular coordinate measured from stag- 
nation, is 

V 

^r=2 sin 6 

For small values of 0 this becomes 

vr 2B -r-m\c) 

Using this value in equation (4), gives 
2.65c 8 ff\ \cj [o 

■ *• 

so that for small values we obtain the approximate value 

_0 I 28cf r 
B e c 

A better approximation is to consider the nose of the airfoil 
elliptical in form. By such an approximation, if the stagna- 
tion point is located at the leading edge, the calculated value 
of 5^ lies between the value obtained for the cylindrical 
nose- and half this value, depending on the thickness ratio. 
Neither the order of approximation nor the importance of 
the value of S^t* warrants much refinement of this calcula- 
tion. It is recommended that the approximation 



c 2 7 
"KR C c 



(6) 



be used for two-dimensional bodies, where r is the radius of 
curvature at the stagnation point which is not necessarily the 
leading-edge radius. 
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For bodies of revolution, by a similar treatment it may be 
shown that equation (6) is, with R c replaced by B L , a fair 
approximation if r is the least radius of curvature at the 
stagnation point. 

As has been noted, the equations (4) and (5) were devised 
on the assumption that the velocity profile of the boundary 
layer was of the Blasius type. In the region behind the 
minimum pressure point the velocity profile changes to a 
separated profile so that if these equations were used to 
determine 5 in this region, it would be expected that the 
calculated would exceed the actual boundary-layer thickness. 
A comparison of the experimentally determined boundary- 
layer thickness for the NACA 0012 airfoil (reference 6) with 
that calculated by this equation and by the more exact but 
laborious method of reference (7) shows (see fig. 1) that this 
equation appears to yield satisfactory results over the 
separating region up to the separation point. 
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Distance along- airfoil surface measured from stagnation point, 
percent chord, s/c x 100 

Fiqtoe 1.— The calculated and measured boundary-layer thickness on the surface of an 
NAOA 0012 aHou at lero lift. 

In the preceding theory the effect of the temperature dis- 
tribution within the heated boundary layer on the boundary- 
layer thickness and velocity profile has been neglected. 
This effect is known to be small for the temperature differ- 
ences necessary for thermal ice-prevention and need not be 
considered for such practical applications. 

DISCUSSION AND CONCLUSIONS 

Recently an experimental investigation of a heat de-icer 
installation on a Lockheed 12-A airplane wing in flight was 
conducted. In the course of this investigation the average 



heat-transfer coefficient from the forward 10 percent of the 
airfoil surface to the air was determined. 

For purposes of comparison the chordwise distribution of 
the heat-transfer cofficient was calculated for the Lockheed 
wing at three spanwise stations and the average heat-transfer 
coefficient determined. The calculated value of 13. Btu per 
hour per square foot per degree Fahrenheit agrees well with 
the value of 11 Btu per hour per square foot per degree 
Fahrenheit as determined from flight tests. 

In the past, for lack of appropriate experimental data, 
heat-transfer coefficients have been calculated on the assump- 
tion that the experimental data of reference 1 for the Clark 
Y airfoil at a given angle of attack may be considered to apply 
to airfoils of other shapes provided the results are corrected 
for the effect of scale found to apply over the limited range 
of test Reynolds numbers for the Clark T airfoil. On this 
basis, the heat-transfer coefficient for the Lockheed 12-A 
airplane wing would be 18 Btu per hour per square foot per 
degree Fahrenheit. The agreement, such as it is, must be 
considered fortuitously close. The assumption that, at the 
same angle of attack and Reynolds number, the heat- 
transfer coefficient for the Lockheed wing (NACA 230-series 
sections) and the Clark T will be identical is clearly unjusti- 
fied. Moreover, there is no basis for the assumption that 
the heat-transfer characteristics at high Reynolds numbers 
may be determined by extrapolation of low Reynolds 
number test results using an extrapolation formula of the 
kind 

o 

where n and a are constants. 

Since the heat-transfer coefficient is detennined by the 
boundary-layer thickness which in turn is dependent on the 
pressure distributioa, it is essential that the pressure distri- 
bution used in the determination of the heat-transfer coeffi- 
cient by the method of this report, must be that corresponding 
to the proper wing Reynolds number. In the event that the 
required experimental pressure distribution is not available, 
the use of the methods of references 8 and 9 combined or the 
more laborious but more exact method of reference 10 as 
modified by the method of reference 11 is recommended. 
Without this modification, the method of reference 10 is 
unsatisfactory at high angles of attack as is indicated in 
reference 11. 

In regard to the application of the calculated wing surface- 
to-air heat-transfer coefficient given by equation (3) to the 
determination of the rate of heat transfer from the airfoil 
surface, it should be noted that the temperature difference 
used in this calculation should be the surface temperature 
minus the air temperature as it is increased as a result of 
fluid friction. The air temperature rise due to fluid friction 

is approximately 1.7 (jM/ ™ degrees Fahrenheit, where V 0 
is the airplane airspeed in miles per hour so that the rate 
of heat transfer is 



(7) 
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where 

Q rate of heat transfer 

S heated area 

T, surface temperature 

To ambient air temperature 

and all variables are in consistent units. The correction 
due to aerodynamic heating is" only important at very high 
speeds and so normally may be neglected. 



Ames Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 

APPENDIX 

Sample calculation: 

In order to illustrate the method of determining the chord- 
wise distribution of the heat-transfer coefficient developed in 
this report, calculations were made for the wing of the Lock- 
heed 12-A airplane used in the experimental flight investiga- 
tion of heat de-icing The wing section used for this illustra- 
tion was taken at a span wise station 123 inches from the 
center line of the fuselage. 

For this calculation the following were assumed: 

Air temperature (° F.) ___ 30 

Altitude (feet) - — 8,125 

True airspeed (miles per hour) . 173 

The wing-lift coefficient was 0.424. At spanwise station 
123 the chord was 7.78 feet and the airfoil thickness was 14% 
percent of the chord. Reynolds number, based upon the 
chord length, was 

i? e -=-^= 10,330,000 

where 

c chord (7.78 ft.) 

V 0 velocity of ambient stream (253 ft. /sec.) 
v kinematic viscosity (0.0001906 ft. s /sec.) 
It was first necessary to calculate the chordwise pressure 
distribution because no data could be found which were 
applicable to this case. The section-lift coefficient was com- 
puted by the method given in reference 12. The chordwise 
normal-force distribution was then found, using the method 
given in reference 8. To find the chordwise pressure distri- 
bution on the upper and lower surfaces of the airfoil, and from 
this the corresponding velocity distribution, the method of 
reference 9 was used. These calculated values have been 

tabulated in the form (1— P) and^pv where V 0 is the veloc- 
ity in ambient stream. (See table I.) A plot of the pressure 
distribution was made to determine the stagnation point 
(fig. 2). The position of the laminar-separation point on the 
upper surface was found by the method of reference 5. 
(See fig. 3.) 
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Figure 2.— The chordwise pressure distribution at wing station 123 of tho Lockheed 12A 
airplane at wing Cr.-0.424 and 10,330,000. 

The heat-transfer coefficient h was calculated by equation 
(3) of this report 

h= ~S 

S, the laminar boundary-layer thickness was computed by 
equation (4) of this report 



. 2 5.3c 3 / 8i/c \ 

5 = ^km~o) 



Substituting value of c=7.78 feet and B c = 10,330,000 
this becomes 



••=0.000081 (^ 



. (W"Ci)\ 



The term 



was evaluated graphically, and is given in Table I. 

Equation (6) was used to find 5 near the stagnation point. 
In the calculation of h the value of X used was 0.765, except 
between the minimum pressure point and the laminar- 
separation point on the upper surface of the airfoil. In this 
region X was reduced linearly from 0.765 at the minimum 
pressure point to 0 at the laminar-separation point. This 
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Fioube 3.— Velocity boundary layer thickness and beat transfer coefficient distribution along the surface at wing station 123 of the Lockheed 12A airplane at wing Cl- 0.424 and R>- 

10,330,000. 



variation of X has been discussed in the section of the report 
giving the .theoretical development of h. 
For these calculations 



1.24X10- J X 



Substituting for k = thermal conductivity of the air at 
30° F. = 1.24 X 10 _s Btu per hour per square foot per degree 
Fahrenheit per foot. 

Values of X, 5, and h are given in table I. 

The chordwise distribution of 5 and h have been plotted in 
figure 3. 
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APPENDIX 

TABLE I. — CALCULATIONS OF THE HEAT-TRANSFER COEFFICIENT FOR WING STATION 123, LOCKHEED 12-A AIRPLANE 



Upper surface of airfoil 


Station 
percent 
chord 


Station 
percent 
chord 


(1-Pri 


V 
V, 


ti 

0.000031-^X10' 


far-© 


to 




ho 


i/eXlOO 


'-xioa 




vi 




(ft.X10») 




(Btt^hr. ft.» 


0.6 
0 

1.26 
2.6 
6.0 
7.5 

10 

16 

20 

30 


0 

.05 
4.75 
6.30 
9.25 
11.90 
14.50 
19.80 
24.70 
34.80 


1.447 

L879 

2.16 

2.21 

2.20 

2.07 

L942 

L743 


L202 
L371 
L468 
1.486 
1.483 
L438 
1.393 
1.322 


i"l88 

1.432 

1.9S8 

2.49 

3.04 

4.23 

5.60 


.244 
.359 
.465 
.569 
.841 
L103 


0.0605 

"."546 
.691 
.839 
1.076 
1.316 
1.888 
2.46 


a765 
.765 
.765 
.765 
.766 
.766 
.651 
.425 
.202 


167 

17.40 
16.08 
11.31 
8.82 
6.14 
2.80 
1.02 


Lower surface of airfoil 


Station 
percent 
chord 


Station 
percent 
chord 


(1-Pi) 


V 
Vo 


i! 

0.000031-^X10» 


f (to"'® 


Sl 


X 


hi 


i/eXlOO 


'-X100 




T, 








(Bta^r. ft.« 


0 

1.26 
2.5 
5.0 
7.5 

10 

15 

20 

30 


"6."9S 
2.35 
5.0 
7.6 

iai 

15. 15 
20.20 

so. 30 


OS 

.752 
.889 
.926 
.962 
1.047 
L112 
1.168 


Oil 
.867 
.943 
.962 
.981 
1.023 
1.054 
1.077 


Oil 

.844 
L640 
2.43 
3.19 
4.65 
6.95 
8.76 


aioo 

.423 
.507 
.601 
.620 
.681 
.562 
.661 


0.451 
.598 
.912 
1.208 
1.407 
1.627 
1.829 
2.41 


0 765 
.765 
.765 
.765 
.786 
.765 
.766 
.766 
.766 


21.1 

15.88 . 
10.42 

7.86 

6.75 

5.84 

6.19 

3.94 
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